Unusual magnetism of layered chromium sulfides MCrS2 (M = Li, Na, K, Ag, and Au) 
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MCrS2 compounds [M — Li, Na, K, Cu, Ag, and Au) with triangular Cr layers show large variety 
of magnetic ground states ranging from 120° antiferromagnetic order of Cr spins in LiCrS2 to double 
stripes in AgCrS2, helimagnetic order in NaCrS2, and, finally, ferromagnetic Cr layers in KCrS2. 
On the base of ab-initio band structure calculations and an analysis of various contributions to 
exchange interactions between Cr spins we explain this tendency as originating from a competition 
between antiferromagnetic direct nearest-neighbor d-d exchange and ferromagnetic superexchange 
via Sp states which leads to the change of the sign of the nearest neighbor interaction depending 
on the radius of a M ion. It is shown that other important interactions are the third-neighbor 
interaction in a layer and interlayer exchange. We suggest that strong magneto-elastic coupling is 
most probably responsible for multiferroic properties of at least one material of this family, namely, 
AgCrS2. 

PACS numbers: 71.20.Lp, 71.70.Gm, 75.30 Et 
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I. INTRODUCTION 

Frustrated magnetic systems attract now considerable 
attentionii Among them there are systems with very 
strong geometric frustrations (e.g. kagome or pyrochlore 
systems), and also less frustrated ones - e.g. systems 
with triangular lattices. Triangular magnets are over- 
constrained and most often they display one or the other 
type of magnetic ordering. Nevertheless, frustrated na- 
ture of triangular layers strongly influences their mag- 
netic properties, often making them rather unusual and 
very sensitive to small variations of the electronic and lat- 
tice structured Such materials also present definite prac- 
tical interest, e.g., as possible thermopower materials^ or 
new multiferroics4i^ 

The presence of orbital degeneracy may introduce spe- 
cial features in the properties of triangular magnets, see 
e.g., Ref. 0. But even without such degeneracy, as in ma- 
terials containing half- filled (i-(sub)shells (Fe^"*" (P; Cr^^ 
tjg), the properties of such systems can be rather non- 
trivial. 

Contrary to similar materials with oxygen instead of 
sulfur, MCrS2 compounds are much less studied. But 
it was recently shown that at least some of them, such 
as AgCrS2, show very interesting behavior: this partic- 
ular material belongs to a pyroelectric class, below the 
Neel temperature Tn = 50 K it develops a rather unusual 
double-stripe (DS) magnetic order— and also becomes 
multiferroic 1^ Motivated by this findings, and trying to 
understand the reasons for this unusual type of magnetic 
ordering, apparently also relevant for the appearance of 
ferroelectricity, we undertook a study of this and similar 
systems with the Af-ions Li, Na, K, Cu, Ag, Au. These 
systems, though in principle very similar and all contain- 
ing as the main building block the same CrS2 triangular 
layers, show very different magnetic ordering: from the 



pure nearest neighbor antiferromagnetism (120° struc- 
ture) for LiCr02^ with the smallest M-cation Li"*" and 
up to ferromagnetic (FM) CrS2 layers in KCrS2-" with 
the largest Af-ion K"*", with more complicated magnetic 
structures in the other systems. Our ab-initio and model 
calculations allow us to explain the general tendency of 
magnetic ordering in this very rich class of compounds, 
and this understanding may be helpful not only for these 
compounds, but also for other magnetic systems with tri- 
angular layers. 



II. CRYSTAL AND MAGNETIC STRUCTURE 

The crystal structure of MCrS2 series has been de- 
termined in Refs. a EMI 

Cr atoms form a triangular 
lattice within CrS2 layers, and the latter are joined by 
M atoms [Fig. [Tfa)]. Cr atoms are located at the cen- 
ter of trigonally distorted octahedra composed of sulfur 
ones. Each S atom is shared by three different octahe- 
dra. But the "connection" between layers is different in 
different compounds. In compounds with alkali metals, 
Li, Na, and K are also sitting in Sg octahedra. One can 
visualize the structure of these compounds as originating 
from the rock-salt structure of (actually hexagonal) CrS, 
in which Cr and alkali ions are ordered in consecutive 
(1,1,1) planes, so that Cr as well as Li, Na, or K are oc- 
tahedrally coordinated by anions (the detailed stacking 
of Cr, S and alkali layers maybe different). 

At the same time, the structure of the systems A/CrS2 
with A/=Cu, Ag, Au is different. In corresponding oxides 
the nonmagnetic ions Cu^^ and Au^^ with the configu- 
ration are linearly coordinated. They are located in 
the centers of oxygen dumbbells, i.e., are sandwiched be- 
tween two oxygens belonging to different MO2) layers. 
The resulting structure is that of delafossitesJ^ 
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FIG. 1. (Color online) High temperature rhombohedral crys- 
tal structures of KCrS2 (a), AuCrS2 (b), and AgCrS2 (c). 
Also shown is a distorted CrSe octahedron surrounding a Cr 
ion in the cell center. 



The structure of their sulfur analogues is more interest- 
ing: it is "in between" that of, say, LiCr02 and AgCr02. 
Such ions are located on top of a S3 triangle of one, 
say, the lower CrS2 layer, but are connected by vertical 
bonds to one S^~ ion of the next, upper layer [Fig.[ljc)]. 
In effect Cu, Ag, and Au are in a "tripod" made of 
four S ions, or in the distorted (elongated in c-direction) 
S4 tetrahedron. Metal ions in such S4 tetrahedra are 
strongly shifted towards the upper, apical S ion. All such 
tripods, or tetrahedra, are pointing in the same direction, 
e.g., up, so that the resulting structure does not have an 
inversion symmetry and is of a pyroelectric class. How- 
ever this interesting structural feature, though probably 
important for some properties of these materials, seem 
to play minor role in magnetic properties of these com- 
pounds, which mainly depend on interactions in CrS2 
layers. Whereas most structural studies of MCrS2 with 
M =Cu and Ag give this structure with M-ions in sul- 
fur "tripods" and R3m symmetry,^ there are also reports 
of a different crystal structure. Thus, in the recent pa- 
per Ref. it is concluded that the symmetry of AuCrS2 
is i?3m or, maybe, i?3m, and the actual structure is the 
delafossite one with linearly coordinated Au+ [Fig.[TJb)]. 

MCrS2 compounds have diverse magnetic structures 
and a broad set of physical properties. Being coupled 
antiferromagnetically (AFM) between the layers, they 
exhibit quite different in-plane ordering at low tempera- 
tures. 

At high temperatures LiCrS2 belongs to P3ml space 
symmetry group. According to neutron scattering mea- 
surements the magnetic structure of this compound be- 
low the Neel temperature Tjv = 55 K exhibits a tri- 
angular spin arrangement (120° structure) within the 
triangular planes, with adjacent planes being coupled 
antiferromagneticallyji^ii^ This structure is typical for 
Heisenberg antiferromagnets with nearest neighbor cou- 



pling on a triangular lattice. The observed value of Cr'^+ 
spin magnetic moment equals 2.26/xb, being considerably 
smaller than the expected value of 3/xb. The difference 
may be presumably attributed to covalency effects, which 
can considerably alter the distribution of the spin den- 
sity around the Cr^+ ion. Indeed, one can expect such 
behavior, keeping in mind much smaller size of Li+ ions 
and respective reduction of unit cell volume. 

KCrS2 undergoes AFM transition at Tjv = 38 KJ^ the 
symmetry group at high temperature is rhombohedral 
R3m. The magnetic structure, in contrast to LiCrS2, 
consists of ferromagnetic layers perpendicular to the c 
axis, which are antiferromagnetically coupled to adjacent 
ones. The paramagnetic Curie temperature of KCrS2 
is not low (0c' = 112K) and indicates that the ferro- 
magnetic interaction in the planes is the dominant one. 
The observed value of the Cr"^"*" spin magnetic moment 
(3.04±0.05/xb) obtained by neutron scatteringi^ is in 
good agreement with the expected value of 3/LtB and with 
the value obtained from the susceptibility measurements 
(3.1/iB). This can be interpreted as an indication that in 
KCrS2 covalency effects are relatively weak. 

In contrast to LiCrS2 and KCrS2, AgCrS2 undergoes 
at Tjv = 41.6 K a first-order phase transition from a para- 
magnetic rhombohedral i?3m structure to an antiferro- 
magnetic monoclinic Cm structurciii Most interesting, 
the material was found to be ferroelectric below T^r, i.e., 
it is a multiferroic system.^ Note that this phenomenon 
differs from the eventual polarization of AgCrS2 due to 
its pyroelectric crystal structure: this polarization ap- 
pears only in a magnetically-ordered state and lies in 
the a6-plane, not along c-direction, as the eventual py- 
roelectric polarization due to the crystal structure it- 
self. In addition to being ferroelectric below Tjv, the 
low-temperature phase of AgCrS2 exhibits an unconven- 
tional coUinear magnetic structure that can be described 
as double ferromagnetic stripes coupled antiferromagnet- 
ically, with the magnetic moment of Cr^^ aligned along 
the h axis within the anisotropic triangular plane. Fer- 
roelectricity below T/v in AgCrS2 is explained as a con- 
sequence of atomic displacements at the magnetoelasti- 
cally induced structural distortion, most probably driven 
by the double-stripe magnetic structure itself. Thus, this 
system can be classified as a type-II multiferroic i^i^^ii 

Similarly to AgCrS2, AuCrS2 undergoes a first-order 
magnetic and structural phase transition at = 47 K 
from a paramagnetic rhombohedral Rim to a monoclinic 
antiferromagnetic C2/m structurej^ The simultaneous 
observation of magnetic and structural transition both in 
AgCrS2 and AuCrS2 gives evidence of a large magnetoe- 
lastic coupling in these systems. This coupling accounts 
for the stability of the observed magnetic order, consider- 
ing that the structural distortions at the transition sup- 
press the geometric frustration of the Cr layers. As we 
will show below, the peculiar antiferromagnetic structure 
observed both in AgCrS2 and AuCrS2 is explained by the 
interplay of the exchange due to direct dd hopping and 
that via anions (sulfur) involving nearest neighbor and 
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further neighbor Cr-Cr interactions, as well as the resid- 
ual frustration in the triangular Cr planes. 

In Table U we put different compounds in the order of 
increasing Cr-Cr distance, which also corresponds to an 
increase of a Cr-S-Cr bond angle since average Cr-S dis- 
tances vary much less than the Cr-Cr ones. One imme- 
diately notices a definite correlation between the crystal 
structure and magnetic order: with increasing Cr-Cr dis- 
tance and Cr-S-Cr angle the magnetic structure changes 
from the 120° AFM structure in LiCrS2 with the small- 
est Li"*" ion and the shortest Cr-Cr distance to ferromag- 
netic layers in KCrS2 with the largest K"*" ion and the 
longest Cr-Cr distance. The crossover between these lim- 
iting cases occurs via incommensurate magnetic phases 
in CuCrS2 and NaCrS2 and the double-stripe structure 
in AuCrS2 and AgCrS2. It is this correlation between 
crystal and magnetic structure, which is the main topic 
of our study. We approach this problem by performing 
ab initio calculations, in which we obtain the electronic 
structure of the MCrS2 compounds, as well as the val- 
ues of relevant exchange constants. We then analyze the 
observed general trends in a superexchange model, dis- 
cussing different relevant, often competing contributions 
to the total exchange. 

TABLE I. Cr-Cr (dcr-Cr) and Cr-S (dcr-s) interatomic dis- 
tances (in A) as well as Cr-S-Cr bond angles 9 (in degrees) 
for the high-temperature A/CrS2 structures. For M=Cu and 
Ag only averaged dcr-s is shown, while 9 is given for two in- 
equivalent S ions. 



M 


rfCr-Cr 


rfCr-S 


9 


magnetic structure 


Li 


3.4515 


2.4063 


91.7 


AFM 120° 


Cu 


3.4728 


2.4036 


90.6, 94.6 


spiral ordering 


Au 


3.4826 


2.3862 


93.7 


AFM double stripes 


Ag 


3.4979 


2.4085 


92.2, 94.1 


AFM double stripes 


Na 


3.5561 


2.4249 


94.3 


spiral ordering 


K 


3.6010 


2.4123 


96.6 


FM in plane 



III. COMPUTATION DETAILS 

Band structure calculations were performed using the 
linear mufBn-tin orbitals (LMTO) method^^ as imple- 
mented in the PY-LMTO computer codc.^*^ We used 
the Perdew-Wang^*' parameterization for the exchange- 
correlation potential in the local spin-density approxima- 
tion (LSDA) . Brillouin zone integrations were performed 
using the improved tetrahedron method. 

When the spin-orbit coupling is not taken into account, 
the use of the generalized Bloch theorem^^ makes possi- 
ble self-consistent calculations of the band structure and 
the total energy E{(i) for spin-spiral structures with an 
arbitrary wave vector q as described in details in Refs.[23l 
andU^- In these calculations the magnetization direction 
in an atomic sphere centered at t + R, where t specifies 
its position in a unit cell and R is a translation vector. 



is defined by two polar angles 9 and </) = q • R -I- (/)o . In 
the present work we considered only planar spin spirals 
with all 9 = 7r/2. The phase (pcr inside spheres surround- 
ing Cr ions was fixed by requiring 0cr = q • tcr, whereas 
for all other spheres it was determined selfconsistently by 
diagonalizing the corresponding spin-density matrix. 

This general approach allows us to treat on the same 
footing not only coUinear, e.g., ferromagnetic or stripe, or 
non-collinear, e.g., 120° AFM, commensurate magnetic 
structures, but also perform calculations for incommen- 
surate helical structures. The only restriction is that it 
should be possible to describe the magnetic structure by 
a single wave vector q. After the q dependence of the 
total energy has been calculated, effective exchange in- 
teractions between Cr spins can be obtained by mapping 
E{q) onto a relevant Heisenberg-like model. 

The magneto-crystalline anisotropy was estimated by 
using the force theorem,-^ i.e., by comparing band ener- 
gies obtained for selected collinear spin structures from 
spin-polarized relativistic calculations with the magneti- 
zation parallel to different crystallographic axes. Spin- 
orbit coupling in these calculations was included into the 
LMTO Hamiltonian at the variational step,^^ 

In order to study the effect of relatively strong elec- 
tronic correlations in the Cr d shell on the band structure 
and magnetic interaction in the AfCrS2 compounds, for 
some of them we also calculated E{q) using the rotation- 
ally invariant LSDA+C/ methodi^ For the double count- 
ing term the so-called atomic limit was used.^^ Other 
details on the implementation of the LSDA-I-C/ method 
in the PY-LMTO code are given in Ref. IH. Calcula- 
tions were performed for the Hund's exchange coupling 
parameter J = 0.9 eV and the on-site Coulomb repulsion 
U = 1.9, 2.9, and 3.9 eV, which gives 1, 2, and 3 eV for 
U,s^U- J. 



IV. RESULTS AND DISCUSSION 

A. Band structure and energies of different 
magnetic structures 

Our band-structure calculations demonstrate that all 
atoms in MCrS2 compounds exhibit their valences corre- 
sponding to the stoichiometry of the compound, i.e., the 
atomic charges correspond to M+, Cr^+, and S^~. The 
s orbitals of M+ are empty, whereas the p orbital of 
are fully occupied. Since Cr atom is triply ionized, there 
are three d electrons localized on a Cr^+ ion. 

The octahedral crystal field at the Cr site causes the 
d orbitals to split into a triplet t2g {xy, xz, yz) and 
a doublet Cg {3z^ — r^, — y^), with the energy of 
the t2g orbitals being lower than that of the Cg states. 
Since there are three d electrons localized on a Cr site, in 
spin-restricted band structure calculations the t2g states 
are half-filled, whereas the Cg levels are empty. In spin- 
polarized calculations the spin-up t2g states are occupied, 
and spin-down t2g are empty. 
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FIG. 2. (Color online) Partial densities of states in AgCrS2 
with the FM alignment of Cr moments. Energies are given 
relative to the Fermi level Ef- 

The cubic component of the crystal field at the Cr 
site is strong enough for the t2g and Cg orbitals to form 
two non-overlapping sub-bands separated by an energy 
gap of about 0.5 eV. Additional trigonal distortion along 
the c-axis lifts the degeneracy of the t2g levels and splits 
them into a singlet aig and a doublet e'g (ai and e' in the 
compounds with i?3m symmetry, e.g., AgCrS2), which 
are linear combinations of the t2g orbitals. Three Cr'^+ d 
electrons occupy spin-up aig and e'g orbitals. 

The electronic structures and density of states (DOS) 
of compounds in MCrS2 series are similar, so to get de- 
tails specific for current calculations, we consider as an 
example the DOSes obtained for AgCrS2 in ferromag- 
netic spin-polarized LSDA calculations (see Fig. [5]). 

The occupied part of the valence band can be subdi- 
vided into several regions. For all M-ions their valence s 



states are empty and d states (if they exist) are totally 
occupied. These valence s and d states do not contribute 
to the electronic density close to the Fermi-energy Ep- 
In AgCrS2 the Ag4d states appear between — 6eV and 
-1.5eV. 

The occupied S^^ ip states form the broad band with 
the width of 6eV between — 6.5eV and — 0.3eV, being 
strongly hybridized both with Ag 4d and Cr 3d states. As 
will be discussed later, this hybridization between Cr 3c? 
and S 3p states is responsible for superexchange along 
Cr-S-Cr and Cr-S-S-Cr paths. According to our band 
structure calculations these materials are insulating even 
in the ferromagnetic state and even without including 
Hubbard's U. For instance for AgCrS2 the LSDA gives 
the energy gap of 0.55 eV. That is, due to their specific 
electronic structure - half-filled t2g subshell and empty 
eg states of Cr'^+ - they would be band insulators (in 
a magnetically-ordered state). When electronic corre- 
lations are accounted for in LSDA-FC/ calculations, the 
occupied majority-spin t2g states are shifted by [/eff/2 to 
lower energies, whereas the unoccupied minority-spin t2g 
and all eg states move ~ C4ff/2 to higher energies which 
increases the values of the gaps. 

Magnetic properties and the electronic structure of 
A/CrS2 compounds are closely related to the occupancy 
of the Cr3d states, which are spread over wide energy in- 
terval from — 6eV to 4eV and form two non-overlapping 
subbands separated by energy gap. Cr eg and t2g orbitals 
form pda- and pdTr-bonds with sulfur p orbitals, respec- 
tively. The hybridization between occupied Cr spin-up 
t2g and Sp-states at — leV and — 0.5eV is clearly ob- 
served. Being rather small below Ep, the hybridization 
between Cr d and Sp above Ep is larger for eg states and 
is well pronounced for spin-up DOSes. 

Our calculations prove the clearly insulating nature 
of these materials. The exchange splitting Aqx ^ 2 eV is 
prominent for the Cr id bands in the whole MCrS2 series 
where only spin-up Cr aig and e'g orbitals are filled. 

The calculated values of the Cr spin magnetic moment 
are close to 3//b for all compounds in the series. Cal- 
culations for spin spirals showed that the Cr moment 
depends only weakly on the wave vector of a spiral, i.e., 
on the kind of magnetic order. In LiCrS2, for instance, 
the moment varies from 2.74^b for the 120° AFM struc- 
ture to 2.98//B for the FM one. This also confirms the 
localized character of the Cr moments and suggests that 
magnetic interactions between them can be described by 
the Heisenberg model. 

Damay et al. in Ref. T analyzed dynamic correla- 
tions and found a small spin gap at very low energies 
as q — > that has been attributed to the weak mag- 
netic anisotropy; i.e., we conclude that the Cr spins in 
AfCrS2 are relatively isotropic and can be described by 
the Heisenberg model. The localized character of Cr'^"'" 
spin magnetic moments is confirmed in our calculations 
by the fact that Cr spin-up aig and e'g states are fully 
occupied, localized on the Cr'^+ site and separated from 
empty states by an energy gap. 
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FIG. 3. (Color online) Representation of double stripe (blue 
arrows) and 90° (red arrows) magnetic structures within the 
Cr plane. One underlying (Si) and one overlying (S2) sulfur 
layers are shown as well. The low-temperature intraplane 
exchanges between the first {Jix, Jiy), the second {J2x, J2y), 
and the third {J^x, Jay) neighbors are shown as curved lines 
with arrows. In the high-temperature phase Jix = Jiy = Ji , 

J2x = J2y = J2, Jix = Jsy = Js- 



The applicability of the Heisenberg model allows us to 
investigate the wide range of Cr spin moment configu- 
rations within the single approach using the Heisenberg 
Hamiltonian in the form 



(1) 



Everywhere below we work in the orthogonal coordi- 
nates, choosing y-axis along one of the directions between 
Cr-Cr nearest neighbors in the a&-plane, and the x-axis is 
chosen perpendicular to it. i.e., it points from one Cr to 
its second neighbor, see Fig.|3l Thus, in our notation the 
q- vectors of magnetic superstructures are given in these 
coordinates, not in the standard vectors of corresponding 
reciprocal lattices. We measure the in-plane components 
of a q- vector in the units of 27r/a and the out-of-plane 
component in 27r/c. 

In case of an arbitrary wave vector q=: [q^, qy, qz) the 
Heisenberg magnetic energy in these coordinates is 



E{q) = ei(q) + e2(q) + eslq) + e.(q) 



(2) 



where (q) are contribution proportional to the exchange 
coupling constants Ji between i-th Cr neighbors within 
the triangular plane (see Fig. [3]). For the undistorted 
high-temperature (HT) rhombohedral structures 

ei(q) — Ji 2cos{Viqxa/2)coii(qya/2) 

+ cos{qya)], (3) 



e2(q) = J2 cos{qxa\/^) 

-f 2 cos(\/3ga:a/2) cos(3gya/2) 



(4) 



e3(q) = J3 2cos(g2;a-\/3) cos((7j^a) + cos(2gya) . (5) 

An expression for interlayer coupling, ez(q), is partic- 
ularly simple for LiCrS2: 



ez(q) = Jz cos(qzc) . 



(6) 



In other compounds with the abc stacking of Cr layers 
Cr neighbors in adjacent planes sit above and below the 
centers of triangles, i.e., above Si and below S2 positions 
in Fig. [31 and (q) becomes 



ez(q) ^ Jz 2 cos(qzc/3 - q^a/ (2^3)) cos(gj,a/2) 



- cos(qzc/3 + q^aj \fZ) 



(7) 



When Jz is sufficiently strong it may affect in-plane mag- 
netic order. 

As sketched in Fig. |31 in the monoclinic low- 
temperature (LT) phases of AgCrS2 and AuCrS2 ex- 
change interactions J„x and J„y between n-th neighbors 
along X and y directions are no longer equal and the 
expressions ([3])" ([5]) should be modified accordingly. For 
instance, the energy of the nearest neighbor interaction 
becomes 



ei(q) = 2JixCO'a{qxX\x)co&{qyy\x) 
+Jiy cos{qyyiy), 



(8) 



where a vector ri^/y (xixiy^yi^/yfi) connects a Cr site 
with its nearest neighbors along x and y directions. 

In order to estimate the effective exchange parameters 
Ji and Jz we first carried out ab-initio calculations for a 
number of q- vectors lying in qz=Q [Fig.HlJa)] and qz='il2 
[Fig. HJb)] planes. The latter value of qz results in 180° 
rotation of Cr spins in adjacent layers. We then fitted the 
q-dependence of the calculated total energy i?(q) (open 
black circles in Fig. S]) by the Heisenberg model given 
by dH)-© using a least-square fit with four (Ji, J2, J3, 
and Jz) and seven {J\x,\y, Jix^iy, J^x.iy, Jz) exchange pa- 
rameters for the HT and LT phases, respectively. The re- 
sults of such a fit for the most interesting system AgCrS2, 
which has the unusual double-stripe magnetic structure 
and becomes multiferroic below T/v, are shown in Fig. |4] 
by filled red circles. A good agreement between the re- 
sults of the LSDA total energy calculation and of the fit 
proves the possibility to describe the magnetic properties 
of these compounds by the Heisenberg model which in- 
cludes the exchange coupling constants between first, sec- 
ond, and third neighbors, plus interlayer exchange con- 
stant Jz- From these calculations we can extract the 
values of the exchange constants for different materials, 
and by comparing the energies of different states we can 
determine which state would be the ground state for one 
or the other system. 



6 



(a) 




(0, 1, 0) 



(2/V3, 0, 3/2) 



q^ff <°'°'^/2) III 



(0, 1, 3/2) 



(c) 



0-7< 



-5 

I -10 

^ -15 
c 

LU 

-20 
-25 



-30-- 







■ \ v '< 

: \\ i Pi / / 


\ \-..-.'-. 

, I-'-'' / 
l"\ i , 




\ ■ ' /} 
\ • ' // 
\ ■ '7 
\ ■ '' f 
\ '■ // 


1 \J\\ 


III V 



(2/V3, 0, 3/2) 



/ (0, 0, 3/2) 



IV (0. 1:3/2) 



FIG. 4. (Color online) Calculated (open black circles) and fit- 
ted (full red circles), using the least-squares method in Heisen- 
berg model, profiles -E(q) of magnetic energies for AgCrS2 
for two sets of wave vectors, see the text. Partial contribu- 
tions in magnetic energy, calculated according to ©-([T}, are 
shown as well. The ei(q), e2(q), e3(q), and ez(q) profiles are 
presented by solid blue, dashed green, dashed with one dot 
red, and dashed with two dots magenta lines. The dispersion 
curves shown in panels (a) and (b) are calculated for the high- 
temperature phase with ferro (a) and antiferro (b) interlayer 
ordering, whereas those in panel (c) are calculated for the 
low-temperature structure and antiferro interlayer ordering. 
The 90° structure, by which we model the double-stripe mag- 
netic structure observed in AgCrS2 and AuCrS2, is marked 
by arrow. 



B. Magnetic structures 

Let us discuss the nature of different competing states 
i.e., different minima of E{q) in Fig. g] The 120° AFM 
spin structure observed in LiCrS2 (Ref. Hsl ) is realized 
by a spin spiral with q= (0,2/3) (in 27r/a units) which 
gives the minimum IV in Fig. |4l The local minimum III 
corresponds to 120° AFM order in the sublattice of 3-rd 
Cr neighbors. 

The FM in-plane structure observed in KCrS2 in 
Ref. [13 corresponds to q = 0. In contrast, the mag- 
netic energy of AgCrS2 (Fig. |4]) has a maximum instead 
of a minimum at this q which agrees with the fact that 
for this system the dominant exchange interactions are 
antiferromagnetic . 

The double-stripe spin structure observed in AgCrS2 
cannot be represented as a single q spiral if the rhonibo- 
hedral unit cell of the HT i?3m structure is used. How- 
ever, it can be easily verified that the Heisenberg energy 
of the DS structure is exactly equal to the energy of a 
spin spiral with q9oo = ('\/3/6,0), shown by red arrows in 
Fig. El in which spins of each i-th Cr chain running along 
the y direction turn by 90° with respect to the previous 
(i — 1) one. 

Indeed, let us consider the interaction of a Cr spin 
from some (j=0) chain with the rest of the Cr plane. 
In the DS structure the spin directions in odd chains to 
the left {-2\i\ + 1) and to the right {2\i\ + 1) are op- 
posite and their contributions to the magnetic energy 
JSq • S2|i|+i=— JSq • S_2|i|+i cancel each other. In the 
90° structure odd chains do not contribute to the mag- 
netic energy because of the orthogonality of Cr spins in 
odd and even chains (JSq • 82^+1 = 0). Consequently, 
the magnetic energy is determined by the interaction of 
So with S2i from the even chains which are exactly the 
same in both spin structures. Here we assume that the 
exchange coupling constants of So with spins from chains 
to the left S_|i| and to the right S|,j| are equal. The cou- 
plings between n-th neighbor lying in the same (Jny) and 
different (Jnx) chains need not to be equal so that the 
degeneracy of the DS and 90° structures holds also for 
distorted Cr layers of the monoclinic LT phase of AgCrS2. 

LSDA supercell calculations performed for the DS and 
90° structures also gave the total energies which are equal 
within the numerical accuracy; with their energy differ- 
ence being less than 1 meV per Cr ion. Because of the de- 
generacy of the two spin structures the energy of the DS 
structure can be calculated within the same spin-spiral 
approach as the energies of other competing magnetic 
states. The corresponding energy minimum is marked as 
II in Fig. H 

Experimentally, however, these two structures would 
lead to somewhat different features of neutron scattering 
spectra although the positions of magnetic Bragg peaks 
are the same. The authors of concluded that the DS 
structure better fits the experimental data than the 90° 
structure. 

The same 90° structure within a Cr plane is also real- 
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ized at q=(3-\/3/6,0) corresponding to the minimum I in 
Fig. m However, because of the rather strong interlayer 
couphng given by Eq. (O the energies at the minima I and 
II are not equaL FinaUy, the maximum at q=(2-\/3/6,0) 
between these two minima corresponds to single-stripe 
magnetic order, in which FM Cr chains running along y 
are ordered antiferromagnetically. 

Comparing the energies of different states in Fig. HI we 
can make several conclusions. First of all, we see that 
if the ordering between planes would be ferromagnetic 
and without extra lattice distortion, Fig.UJ^a), the abso- 
lute minimum for the parameters calculated for AgCrS2 
would correspond to the simple 120° AF structure, i.e., 
the state IV in Fig. S] Such in-plane ordering is indeed 
realized in LiCrS2, but for real AgCrS2 the observed or- 
dering is different and corresponds to the double stripe 
structure. 

When we change the interlayer ordering, making it an- 
tiferromagnetic, the situation already changes: the 120° 
state (state IV) is destabilized, and another state. III, be- 
comes the absolute minimum. Fig. Sl^b). We also notice 
that AFM interlayer ordering strongly lowers the energy 
of the 90° structure (state II in Fig. HJ^b)) which, as dis- 
cussed above, is degenerate with the DS one, so that this 
state starts to compete with the state III. And when we 
include the lattice distortion present in AgCrS2 in the 
LT phase. Fig. SKc) , the double stripe state II becomes 
the absolute minimum. Thus, we see that for the lat- 
tice corresponding to the real LT structure of AgCrS2 
the double-stripe magnetic ordering with the antiferro- 
magnetic coupling between layers is indeed the ground 
state in our calculations. We also see that several factors 
are important for the stabilization of such DS structure: 
besides particular ratio of different exchange constants, 
see below, also particular 31? interlayer ordering and lat- 
tice distortion, accompanying magnetic ordering, are all 
important for making double-stripe structure. 

But we also see from Fig.|3]that there exist, especially 
in the HT lattice, other magnetic states competing with 
the double-stripe one. Thus, one can predict that the 
magnetic fluctuations above T^r, which could be probed 
e.g. by inelastic neutron scattering, could be most pro- 
nounced not at the wave vector corresponding to the 
double-stripe ground state structure, but at other values 
of q, for instance those corresponding to the solutions III 
and IV in Fig. S) 

Yet one more conclusion which we can extract from 
Fig.Uis that, at least in AgCrS2, the spin-lattice (magne- 
tostriction) coupling is very important in these systems: 
only when we included the lattice distortion, occurring in 
AgCrS2 below T^r, did we obtain the real double-stripe 
structure as a ground state. 

In addition to the nonrelativistic calculations discussed 
above, we have also studied the magneto-crystalline 
anisotropy in MCrS2 by accounting for spin-orbit cou- 
pling in calculations for the FM spin structure with the 
magnetization directed along different crystallographic 
axes. It turns out that Cr atoms form an easy- 



plane magnet, which is consistent with the experimen- 
tal results the spin-orbit coupling rotates all Cr 
spin magnetic moments into the ab-plane even in high- 
temperature phase, but does not affect the magnetoelas- 
tic in-plane coupling and low-temperature lattice distor- 
tion. 



C. Exchange constants 

The LSDA exchange parameters estimated for the HT 
structure of all six AfCrS2 compounds by fitting corre- 
sponding E{q) using the Heisenberg model defined by 
©-([T]) are presented in Table HT] We first do not con- 
sider the LT phases of AgCrS2 and AuCrS2, because 
we want to concentrate on general trends observed in 
this whole class of materials. Detailed results for the LT 
phases will be presented below. The dependence of the 
exchange constants on U in LSDA-I-C/ calculations will 
be discussed in Sec. IIV El 

From the Table HIl we see that, with the exception of 
AuCrS2 which deviates from the general trend and will 
be discussed below, the variation of the nearest-neighbor 
exchange Ji in the series M = Li, Cu, Ag, Na, K clearly 
correlates with the corresponding structural parameters 
from Table HI With the increase of the size of M ion and 
of the Cr-Cr distance, Ji changes form strongly antifer- 
romagnetic in LiCrS2, with the smallest Li and short- 
est dcr-Cn to strongly ferromagnetic in KCrS2, with the 
largest K and longest dcr-Cr, and becomes very small in 
the Cu and Ag compounds with intermediate Cr-Cr dis- 
tances. 

We also notice that in all the compounds the third 
neighbor exchange J3 is antiferromagnetic and rather 
strong. On the other hand, the second neighbor exchange 
J2 is weak and can in most cases be neglected. Appar- 
ently it is an interplay of the nn exchange Ji and the 
third neighbor exchange J3 which is primarily responsi- 
ble for the stabilization of one or the other spin structure 
in the AfCrS2 series. 



TABLE II. Different exchange coupling constants (in meV) in 
the high-temperature phase of MCrS2 calculated in LSDA. 



M 


Ji 


J2 


J3 




J1/J3 


Li 


5.17 


0.46 


2.73 


0.93 


1.9 


Cu 


0.16 


0.03 


1.51 


0.82 


0.1 


Au 


7.41 


1.63 


5.93 


2.93 


1.3 


Ag 


-0.14 


-0.13 


2.45 


0.74 


-0.2 


Na 


-4.06 


0.23 


2.49 


0.09 


-1.6 


K 


-5.45 


0.19 


2.11 


0.05 


-2.6 



Taking these considerations into account it seems rea- 
sonable to apply the J1-J3 model to investigate mag- 
netic ordering in MCrS2. It is well known that the 
simple_Ji model with antiferromagnetic Ji > (see e.g., 
Ref. 30) gives noncollinear magnetic ground states with 
q=(0, 2/3) and angles of 120° between spin magnetic 
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moments. In the J1-J3 model the magnetic energy 
equals £'i^3(q) — ei(q) + e3(q). A simple analysis shows 
that for positive Ji and J3 the wave vector qiv = (0, 
2/3) is still the global minimum with the energy of 
E — — 3/2(Ji + J3). Here we consider only extrema at 
wave vectors lying on x and y axes. Other symmetrically 
equivalent extrema can be obtained by applying ±27r/3 
rotations to corresponding q. The numbering of the min- 
ima corresponds to the notations in Fig. |4l 

For I Ji| <4J3 a local minimum appears at qn — {Qx, 0) 
on the X axis, with defined by cos(\/37rga;) = — J1/4J3. 
When Ji < J3/2 another minimum qm appears also on 
the y axis which becomes the global minimum for FM 
Ji < 0. If Ji = 0, qm = (0,1/3) corresponds to 120° or- 
der of 3-rd neighbor spins. As the strength of FM Ji 
increases, both qn and qni shift towards zero, until for 
FM I Ji I > 4 J3 the two minima merge at q = which be- 
comes the global minimum. 

These additional minima at incommensurate qn and 
qm imply possible formation of helical magnetic order, 
but the exact picture does depend on interlayer exchange 
coupling Jz too,.-'^'* the latter being one of possible way to 
stabilize the magnetic structures observed in the "inter- 
mediate" systems MCrS2 (M=Cu, Au, Ag, Na). In par- 
ticular, this may be the origin of incommensurate mag- 
netic structures for M =Cu, Na, or commensurate double 
stripes for M=Ag, Au. 

An extra complication is introduced by the observed 
monoclinic distortion in AgCrS2 and AuCrS2, which in- 
duces three pairs of nonequivalent nearest-neighbor ex- 
change couplings {Jix, Jiy), {J2x, -hy), and [J^x, J^y) 
(see Fig. [3]). The observed four-sublattice spin arrange- 
ment cancels the effect of Jix and J2y In order to clarify 
which of remaining magnetic exchanges are relevant for 
the stabilization of the DS structure, namely the ferro- 
magnetic first-neighbor coupling Jiy, the antiferromag- 
netic second-neighbor J2X and antiferromagnetic third- 
neighbors Jzx and J^y superexchanges, and the inter- 
plane antiferromagnetic superexchange Jz we calculated 
the energy of different Cr spin moment configurations 
and derived the corresponding exchange values. The re- 
sults are summarized in Table IIIII They show that the 
monoclinic distortion does stabilize the DS structure by 
strongly suppressing the AFM contribution to Jiy along 
the FM Cr chains. 



TABLE III. LSDA exchange coupling constants (in meV) for 
low temperature phases of AgCrS2 and AuCrS2. 



M 


Jlx 


Jiy 


J2x 


J2y 




Jay 


J zx 


Jzy 


Au 


5.14 


2.69 


0.70 


0.57 


3.19 


3.24 


1.85 


1.32 


Ag 


1.12 


-1.36 


-0.30 


-0.23 


2.54 


2.62 


1.09 


0.62 



We also have to comment on the values of exchange 
constants for AuCrS2 shown in Tables |ll] and IIIII These 
values definitely deviate from the regularities observed 
in other materials of this series. The ratio of the im- 
portant exchange constants J3 and Ji for AuCrS2 is still 



such that it gives the double-stripe structure observed 
experimentally. However the absolute values of these ex- 
changes for this system are about two times larger than 
what one would expect from the comparison with other 
materials of this class. We do not have a full explanation 
of this difference. A possible reason is that AuCrS2 has a 
delafossite structure with interlayer Au"'' ions in a linear 
coordination.^^ It is possible that the reason for differ- 
ent values of exchange for this system is connected with 
that. Still, this situation is definitely unsatisfactory, and 
it requires further study. 



D. Interpretation of magnetic properties 

Our calculations, presented above, have shown that 
indeed the observed types of magnetic ordering in Cr- 
plane in MCrS2 (120° for Li, double stripes for Ag and 
Au, ferro layers for K) are reproduced. The obtained 
values of exchange constants. Table [Hi allow to explain 
these magnetic structures. 

Thus for the smallest Af-ion Li the nn exchange 
J I is the strongest and antiferromagnetic; apparently 
it is predominantly responsible for the observed pure 
antiferromagnetic (120°) ordering, observed in LiCrS2. 
With the increasing of Cr-Cr distance and Cr-S-Cr 
angle (Li— >Cu— >Au— >Ag->Na->K) the value of Ji de- 
creases and then changes sign, becoming ferromagnetic 
for Af=Ag, K. Simultaneously the AF exchanges between 
3-d neighbors J3 remains relatively large, and it plays 
important role for intermediate compounds Ag/AuCrS2, 
apparently leading to their double-stripe ordering. Fi- 
nally, large nearest neighbor ferro interaction Ji for the 
large M-ion K guarantees ferro ordering in Cr-plane in 
KCrSz. 

To understand the microscopic origin of different ex- 
change integrals in this series, one should look at different 
microscope exchange passes. In Fig. [Sl[a)-[Sl[d) we show 
the main paths of superexchange, existing in CrS2 planes 
with the Cr'^^ ions with d-shell tl^^e^ and with the geome- 
try of edge-sharing CrSg octahedra with nearest neighbor 
Cr-S-Cr angle of about 90°. 

First of all, there exist a direct overlap of different 
t2g orbitals of neighboring Cr ions, e.g., xy-orbital in 
Fig-E^a). It gives a rather large AF exchange 

which however strongly decreases with the increasing of 
Cr-Cr distance. 

In Fig. [S^b) and [SJc) we show an exchange of ^29-^23 
via 90° Cr-S-Cr bond. The processlSfb) (virtual hopping 
of t2g electrons through the same ligand p-orbital, in this 
case pz), gives strong AF exchange: 
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1-st nearest Cr neighbours 
with edge sharing 




Crl\{ Cr2 



(d') ^ (d") 

FIG. 5. (Color online) All possible contributions to Ji 

where we denoted by A the charge-transfer en- 
ergy (the energy of a transition Cr'^+((i'^)S^^(3p^) — > 
Cr^"'"((i'')S~ (3p^)). One sees that this process does not 
change strongly with the Cr-S-Cr angle, only the dis- 
tance Cr-S determines the value of tpd-K hopping, and 
this distance is approximately constant in the whole se- 
ries MCrS2. 

The process [5jc) (the ^29-^29 exchanges via different 
Sp-orbitals) leads to the ferromagnetic exchange, which 
is however usually weaker, 



Jc 



(11) 



(here Jh,s is the Hund's rule coupling on sulfur) and 
it decreases by absolute value with decreasing Cr-S-Cr 
angle. 

More important is another ferromagnetic contribution 
due to a virtual hopping from the occupied i2g-shell of 
one Cr to the empty Cg-shells of another Cr, Fig. [Sjd). 
As one sees from Fig. Eld) this process also gives ferro- 
magnetic contribution. 



4-2 f2 
pd(T pd-TT 

^2 ^2 

pda pdTT 
A3~ 



3Jg 

Udd 

A ' 



AFM 




Cr-Cr 



(12) 



FM 



FIG. 6. (Color online) Schematic dependence of different con- 
tributions to the nearest neighbor Cr-Cr exchange Ji. 



where the first term corresponds to a process [H^d') (ef- 
fective transfer of an electron from one Cr to the other 
via S), and the second EJd") (the transfer of two 3p- 
electrons of S to the left and right Cr ions). We do not 
keep here some numerical coefficients. Note that despite 
the presence of a small factors or this ferromag- 
netic contribution (fT2|) is comparable with ([9|) (typically 
tpda ^ V^tpdTr), and also the Hund's rule contribution 
in (|12l) is enhanced by factor 3. Thus, though usually 
the 90°-exchange involving Hund's rule interaction gives 
ferromagnetic, but weaker exchange, in this case due to 
a specific electrons occupation of Cr^+ it can give sig- 
nificant contribution and can even start to dominate if 
the other competitive contributions are small. This is 
apparently what happens in KCrS2, in which the main 
competing AF exchange [Sfa) is strongly reduced due to 
a large size of K+ and corresponding increase of Cr-Cr 
distance. 

Thus, we can schematically present different contri- 
butions to the nearest neighbor exchange Ji and their 
change in the row (Li— j'Cu— J'Au— ^Ag— >Na— >-K)CrS2 as 
following. Fig. [6l 

To explain the resulting magnetic structures, especially 
double-stripe structure of AuCrS2 and AgCrS2, we also 
have to include the further neighbor exchange. As seen 
from Table HIl the second neighbor exchange is always 
small. Somewhat surprisingly, larger and more important 
turns out to be the interaction of third neighbors. It can 
be schematically explained by the consideration shown in 
Fig. [71 in which one sees that there is an exchange path 
connecting occupied t2g orbitals on third neighbors Cri 
and Cra via two sulfur Si and S2, (with their p-orbitals 
being relatively large) due to the p-p overlap, or to the 
overlap via an empty e^-orbital (a;^ — y^) of Cr2 (Fig. [7]). 

Thus in this geometry the coupling between third 
neighbors J3 turns out to be reasonably large (larger 
than J2) and antiferromagnetic, and in effect it is this 
coupling which stabilizes double-stripe structure for "in- 
termediate" composition AgCrS2 and AuCrS2, in which 
the main nearest neighbor interaction Ji is small due to 



10 




FIG. 7. (Color online) A possible exchange path contributing 
to the antiferromagnetic exchange of third neighbors J3. 

a compensation of different contributions to it. 

The general tendency showing regular change of dif- 
ferent exchange contributions, especially of the nearest 
neighbor exchange Ji, see Fig. [51 is also confirmed by 
model calculations in which we took LiCrS2 and artifi- 
cially compressed it in c-direction, keeping the volume 
constant. At this change the in-plane Cr-Cr distance 
and Cr-S-Cr angle increase, following the same trends as 
in going from LiCrS2 to (Ag, Au) and to KCrS2. Our 
ab-initio calculations of this model system confirmed the 
trend discussed above: with increasing Cr-Cr distance 
large AF coupling Ji strongly decreases and becomes fer- 
romagnetic. 

E. The effect of LSDA+?7 on calculated exchange 
constants 

So far we discussed only exchange coupling constants 
determined by fitting ^^(q) curves calculated within 
LSDA (Table |IT|. Comparing the J1/J3 ratio from Ta- 
ble |lT] with the critical values obtained from the analysis 
of the J1-J3 Heisenberg model one notices that for some 
of the compounds the estimated Ji do not give an experi- 
mentally observed ground state. For KCrS2, for example, 
J1/J3 = —2.6 > —4 corresponds to an incommensurate 
spin-spiral structure in the a-b plane instead of experi- 
mental FM ordering. 

One of possible reasons for this is that the LSDA un- 
derestimates the Coulomb repulsion between rather lo- 
calized Cr 3d electrons. The Udd parameter in expres- 
sions ([^- P^ is the energy cost of adding an electron 
to one of the unoccupied minority-spin t2g state. In 
LSDA it is governed solely by the exchange splitting of 
about 2.4 eV between the minority- and majority-spin t2g 
states, i.e., by Hund's coupling of 3Jh- As a result the 
LSDA overestimates those contributions to the inter-site 
exchanges that have Udd in the denominator. Accounting 
for the Coulomb repulsion in LSDA+U calculations in- 
creases the energy difference between the minority- and 
majority-spin Cr t2g states by Ucs, so that Udd becomes 
equal 3Jh + Ucb- 

The increase of Uctf suppresses AFM Ja and Jb, 
whereas the FM t2g-eg contribution Jd is much less af- 
fected. Thus, in the compounds with FM Ji (M = Na, 



K) it becomes even stronger, whereas in those compounds 
for which LSDA gives AFM Ji its value decreases and it 
may even change sign. On the other hand, the AFM 3- 
rd neighbor coupling J3, which is governed by the t2g-t2g 
superexchange (Sec. lIVDj) . gradually decreases with the 
increase of Ucs- 

This combined effect of strengthening the FM Ji and 
weakening the AFM J3 leads to a reduction in the Ji / J3 
ratio estimated for KCrS2 from —2.6 in LSDA to —4.6 
and —6.4 in LSDA-|-[/ calculations with UcS = 1 and 2 
eV, respectively. Thus, accounting for Coulomb repulsion 
stabilizes the FM in-plane order in KCrS2. In LiCrS2 
the 120° structure gives the lowest total energy also in 
LSDA-|-[/ calculations. In other compounds the increase 
of UcS changes the J1/J3 ratio and, consequently, the 
position of incommensurate minima. 



V. SUMMARY 

Summarizing, the results of our ab-initio calculations, 
and model considerations of Sec. IIVDI allowed us to ex- 
plain the very interesting sequence of magnetic phases in 
layered chromites MCrS2 with triangular Cr layers, in 
which the magnetic ordering in Cr layers changes from 
purely antiferromagnetic (120°) structure in LiCr02 via 
"intermediate" double-stripe structure of AgCrS2 and 
AuCrS2 (and incommensurate structure in NaCrS2 and 
CuCrS2) to a ferromagnetic layers in KCrS2. These 
structures emerge mainly as a result of competing con- 
tributions to the nearest neighbor exchange Ji, together 
with reasonably large antiferromagnetic exchange for 
third neighbors J3. In particular, their combined action 
leads to the most interesting double stripe structure of 
AuCrS2 and AgCrS2, which apparently is responsible for 
the multiferroic behavior of the latter (and probably also 
in the former — it is not checked yet). Our study demon- 
strates quite nontrivial interplay of lattice geometry and 
orbital occupation in giving such diverse magnetic behav- 
ior in apparently rather similar materials. The frustrated 
nature of the lattice definitely plays a very important role 
in these phenomena. Such high sensitivity of magnetic, 
and apparently some other, e.g., multiferroic properties 
to fine details of electronic and lattice structure could 
probably be used also to tune the properties of other 
similar materials. We envisage that further studies of 
the stability of nuclear and magnetic structures may pro- 
vide a clue to tailor the magnetoelastic coupling and the 
multiferroic properties in geometrically frustrated oxides, 
sulfides and selenides with different transition metals. 
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